Frontiers in Neuroimaging

Much of a typical individuals mental life is characterized by spontaneous thoughts that occur independently of external stimuli. In prior studies, ongoing mental experiences and their neural correlates have been captured using thought probes presented at random intervals during functional Magnetic Resonance Imaging (fMRI). However, this approach results in temporally imprecise estimates of brain activity relative to the arising of mental experience. In this preregistered, proof-of-concept study, we aimed to improve temporal precision using a novel method termed real-time fMRI-triggered experience-sampling (rt-fMRI-ES). We analyzed blood-oxygenation-level-dependent signals in real time during a wakeful resting state (n=60) to trigger thought probes from spontaneous activations within two regions: the dorsal anterior insular cortex (daIC; a key region within salience network) and posteromedial cortex (PMC; a key region within default mode network). We tested two preregistered hypotheses: (H1) Ratings of arousal time-locked to daIC-activation trials are higher than ratings time-locked to non-daIC-activation trials; (H2) Ratings of external-attention time-locked to PMC-activation trials are lower than ratings time-locked to non-PMC-activation trials. After applying preregistered exclusion criteria, 42 participants (1243 trials) and 49 participants (1429 trials) were included in H1 and H2 analyses, respectively. We did not find evidence in support of H1, but we did find evidence in support of H2, as external-attention ratings were significantly lower for trials triggered by PMC activation compared to other trial types. Taken together, we successfully developed and validated the rt-fMRI-ES method, offering a novel technique to efficiently capture spontaneous thoughts based on ongoing neural activity. Preregistered Stage 1 Recommendationhttps://osf.io/sd4hu (Date of in-principle acceptance: 07/24/2024; under temporary private embargo)

One of the prominent methods in neuroimaging data processing is SSM-PCA, which is based on principal component analysis and allows for the identification of diagnostically significant patterns in the form of statistical maps. We developed software, PIE Toolbox, employs SSM-PCA and classification based on the obtained diagnostic patterns revealed from functional and structural tomographic brain imaging. The program supports the entire analysis pipeline including preprocessing of brain images, diagnostic patterns extraction, building classification models, and prediction based on them. The resulting diagnostic patterns are weighted principal components obtained through SSM-PCA, or their linear combinations. PIE Toolbox allows selection of relevant structural and functional brain patterns, computation of their expression values in regions of interest, classification using support vector machines, and evaluation of model performance via cross-validation. This approach enables the use of patterns as features of intergroup differences for individual diagnosis. The software has been validated on both simulated and ADNI datasets.

There is growing interest in identifying brain function underlying adolescent cognition, personality, and psychopathology. One promising approach is Precision Functional Mapping (PFM) of MRI functional connectivity, a data-intensive method for characterizing individualized brain networks. Foundational studies suggest that PFM can detect stable, task-responsive, and clinically relevant networks. Studies demonstrate that both functional connectivity reliability and network stability improve with increasing data quantity, although benchmark estimates vary across populations, preprocessing pipelines, and MRI acquisition approaches. Accordingly, it is important to understand how PFM performs in adolescent populations and with multi-echo fMRI acquisition. In a case study of eight youth (ages 10-17), we applied PFM to 80-minutes of combined resting-state and task-based fMRI. The resulting networks were highly modular, consistent with adult templates, and without evidence of structural registration artifacts. Functional connectivity reliability compared favorably to prior single-echo studies, with multivariate similarity and ICC estimates showing early stabilization around 10-15 minutes despite continued improvement with additional data. Trait-like stability increased gradually with acquisition time and a Bayesian algorithm (MS-HBM) demonstrated higher stability than Infomap. Across algorithms, stability was greatest in sensory networks (somatomotor, auditory, visual). Furthermore, when evaluating task-based responses to threat and attention paradigms, only the auditory network consistently benefited from individualized mapping over group template networks. These findings suggest that, with constrained scanning time, PFM is especially effective for characterizing sensory and perceptual networks in adolescents. Bridging the methodological divide between deeply sampled individual cases and large-scale developmental studies will require further innovation and validation.

BackgroundMagnetic resonance imaging (MRI) is critical for acute stroke triage, but time-consuming, and often requires contrast injection for perfusion imaging. This study aimed to synthesize T-map perfusion maps from routinely available, non-contrast DWI and FLAIR using deep generative models. We hypothesized that relevant perfusion information could be inferred from these modalities to streamline imaging and reduce reliance on dynamic susceptibility contrast perfusion. MethodsAcute MRI data from 355 patients with anterior circulation stroke, including dynamic susceptibility contrast perfusion, were retrospectively collected from two European centers (Heidelberg: 2010-2018; Bordeaux: 2021-2022). Six versions of a denoising diffusion probabilistic model (DDPM) and a GAN architecture were trained to generate synthetic T-max perfusion maps from DWI, FLAIR, and infarct core mask as inputs. Performance was assessed by comparing synthetic and ground truth T-max maps using image similarity metrics. Regions with T-max >6s were compared using Dice coefficients, and mismatch volume distributions were analyzed. An ablation study quantified the contribution of each input. ResultsThe best performance was achieved by a DDPM with a 2.5D architecture using DWI, FLAIR, infarct core mask, and a perfusion-weighted loss function. It produced synthetic perfusion T-max maps with high similarity to ground truth under 110 seconds. The model showed strong spatial overlap for T-max >6s regions in internal validation (average Dice = 0.82, SD = 0.08), and external validation average (Dice 0.59, SD = 0.13), respectively. Synthetic maps closely matched ground-truth mismatch distributions, capturing key perfusion patterns. The infarct core mask played a critical role in model performance, alongside DWI and FLAIR inputs. ConclusionsWe propose a non-invasive, scalable framework to generate synthetic T-max perfusion maps from non-contrast MRI. This approach could expand access to perfusion data in acute stroke, shorten imaging protocols, and accelerate treatment decisions by eliminating the need for contrast-enhanced acquisition. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/684079v2_ufig1.gif" ALT="Figure 1"> View larger version (94K): org.highwire.dtl.DTLVardef@164235forg.highwire.dtl.DTLVardef@14e5489org.highwire.dtl.DTLVardef@190214eorg.highwire.dtl.DTLVardef@17a9e3a_HPS_FORMAT_FIGEXP M_FIG C_FIG

IntroductionPostmortem human brain magnetic resonance imaging (MRI) offers a unique opportunity to study finer neuroanatomical details and enables direct correlations with gold standard histological and immunohistochemical assessments. However, to prevent tissue decay, postmortem brains are preserved in fixative solutions which can alter tissue properties and exert substantial impacts on the MRI signals. The present study investigates the impact of formalin fixation, the most commonly used solution for postmortem human brain preservation, on different quantitative MRI contrasts. Methods142 intact human brain hemispheres immersed in 10% formalin for a range of fixation durations (between 0 days and 20 years) were imaged in a 3T MRI scanner. A subset of 10 brains were further scanned repeatedly at days 0, 3, 10, 20, 30, 60, 90, and 120 to allow for better characterization of the initial transient effects of fixation. Voxel-wise T1 and T2* relaxation, T1/T2 ratio, and myelin water fraction (MWF) maps were generated for each specimen and timepoint, and linear and nonlinear models were used to examine the spatiotemporal changes associated with progressive fixation. ResultsAll investigated metrics were significantly impacted by formalin fixation, albeit at different rates and with differing regional patterns. T1 and T2* relaxation time decreased as a result of progressive fixation, whereas T1/T2 ratio and MWF measures increased. T1 relaxation and T1/T2 ratio showed nonlinear patterns with initially accelerated changes that decelerate in the first few months, whereas T2* relaxation and MWF changes followed a more linear trend. ConclusionFormaldehyde fixation exerts systematic changes on quantitative MRI signals that can be modeled and adjusted for to allow for harmonized comparisons of MRI metrics across brains fixed for differing durations. The distinct temporal trajectories observed across metrics highlight the need to account for fixation duration in study design and downstream analyses, particularly when integrating datasets acquired under heterogeneous conditions. Our findings provide a quantitative framework for correcting fixation-induced biases, thereby improving the interpretability and reproducibility of postmortem MRI studies.

Task-evoked decreases in blood-oxygenation-level-dependent (BOLD) signals are a well-recognized phenomenon in functional magnetic resonance imaging (fMRI) studies. These deactivations are most prominent in the default mode network (DMN), a set of regions most active at rest. The metabolic basis of task-induced BOLD fMRI deactivations remains unclear. To address this question, we used PET/MRI to simultaneously measure BOLD fMRI and cerebral glucose consumption (CMRglc) during visuomotor and language tasks in 22 cognitively unimpaired older adults (15 female, 7 male). Task performance increased BOLD signals in task-relevant regions and decreased BOLD signals in the DMN. Positive BOLD responses generally coincided with increases in CMRglc. In contrast, CMRglc did not decrease in regions showing negative BOLD responses; instead, it typically increased. In particular, the posterior cingulate cortex showed significant CMRglc elevations in conjunction with negative BOLD responses. Whole-brain intensity normalization partially restored task-induced decreases in CMRglc, indicating that relative reductions appear in regions in which CMRglc increases are smaller than the global average. Overall, our results imply that BOLD fMRI deactivations can occur in conjunction with stable or even increased glucose consumption.

We trained a self-configuring nnU-Net model for CMB segmentation in a heterogeneous multicenter sample (n=264), including 1.5T and 3T field strengths, SWI and T2*-GRE sequences, and community and clinical cohorts. Model performance was evaluated using 5-fold cross-validation with a focus on object-level detection metrics. Real-world performance was evaluated on scans from an unseen dataset of people with cerebrovascular disease (n=20). The model achieved 0.82 cluster Dice, 0.88 precision, and 0.77 sensitivity on hold-out test data. Notably, the model demonstrated a low false-positive rate, averaging 0.58 false positives (FPs) per scan, an improvement on existing publicly available models. The model achieved high performance in dataset of those with Alzheimer's disease and mild cognitive impairment (0.89 cluster Dice, 0.94 sensitivity), supporting its utility in clinical settings where ARIA-H monitoring is critical. In external validation, the model maintained high robustness with 0.79 sensitivity and 0.95 FPs per scan. By leveraging a heterogenous training strategy and a self-adapting architecture, we demonstrate that deep learning can achieve high-precision CMB detection that is robust to domain shifts. The low FP rate suggests this publicly available pipeline is suitable for automated screening and lesion counting in heterogenous large-scale clinical trials, reducing the burden of manual quantification.

PurposeTo present a first-of-its-kind artificial intelligence (AI-)integrated MR pulse sequence that detects out-of-voxel (OOV) artifacts in real-time (within-TR) and responds prospectively by updating the crusher gradient scheme. MethodsPer Excitation Real-time Execution & Guided Responses with Integrated Neural-network Evaluation (PEREGRINE), developed for deployment of deep learning models and sequence updates, operated time-domain (TD) and frequency-domain (FD) convolutional autoencoders that detect OOV artifacts. Scans without (AI-off) and with (AI-on) updates were collected from the prefrontal cortex of healthy volunteers using edited MRS. The degree of OOV contamination (OOV score) was quantified per transient based upon the prevalence of OOV signals in the TD and FD data. OOV scores above a user-defined threshold triggered an update of the gradient scheme, iterating through 48 permutations (6 axis transpositions x 8 polarity flips). ResultsWithin each 2-second TR, PEREGRINE successfully provided single-transient OOV scores and updated gradients accordingly. No difference was observed between the OOV scores from the full ("Full" condition) AI-on and AI-off sessions due to the AI-on scan cycling over better and worse gradient permutations relative to the AI-off scan. However, the AI-on scan had significantly lower OOV scores than the AI-off scan when selecting the transients where PEREGRINE persisted ("Dwell" condition) on a given gradient permutation. Ultimately, Fit Quality Number (FQN), from linear combination modeling, improved significantly for the AI-on compared to the AI-off scan. ConclusionPEREGRINE enabled an AI-integrated sequence allowing for real-time evaluation and reduction of OOV artifacts, identifying gradient modifications that produced less OOV contamination.

Sensorimotor beta activity (13-30 Hz) is a key neuronal signature in the human sensorimotor system, and its features can be effectively measured using functional brain imaging methods such as magnetoencephalography (MEG). In addition to its importance in healthy brain processing, beta activity has been shown to be altered in several neurological diseases, underscoring its potential as a biomarker. To serve as biomarkers, features must be reliably defined, stable across measurements and, ideally, amenable to automated analysis, yet current approaches to beta characterization require subjective decisions and manual work. We here describe a hidden Markov model (HMM) based approach to automatically segment beta events from source level MEG beta band activity into discrete high- and low-beta states. We demonstrate the differences between the proposed HMM based approach and a commonly used amplitude-envelope based approach to analyse high- and low-beta modulation. We show that the methods complement each other both when applied to resting data and task related passive movement data. Furthermore, we assess the test-retest reliability of the proposed pipeline within individuals using intraclass correlation coefficients (ICC), and test if HMM constructed at one measurement site can be applied to data acquired at another site, thereby evaluating its multisite transferability. We show that the proposed approach produces stable results within subjects and across sites for many of the features. The ICC values were excellent for high-beta state (86-100% of brain areas), while low-beta state test-retest reliability was more modest. Most of the features showed statistically significant differences between sites only in a few brain areas, indicating very good multisite stability. The proposed approach can serve as an automated, reproducible analysis pipeline for, e.g., clinical applications, and appears suitable for multi-site datasets.

IntroductionCoronary artery disease (CAD) increases the risk of cerebrovascular events, yet early brain injury in this population remains poorly characterized. White matter hyperintensities (WMHs), a biomarker of cerebrovascular lesions, are prevalent in CAD and are linked to risk of stroke. Beyond total burden, spatial distribution of WMHs carries pathological significance and is critical for understanding CAD-related injury. While clinical outcomes including coronary revascularization procedure and myocardial infarction influence CAD prognosis, their impact on WMH burden remains unclear. MethodsThis study investigated regional WMH burden in CAD and its relationship with clinical characteristics. 82 adults over 50 years participated, including 44 individuals with CAD and 38 controls. WMHs were segmented from fluid attenuated inversion recovery and T1-weighted MRI and categorized as total, periventricular, deep, and superficial regions. History of myocardial infarction and coronary revascularization (coronary artery bypass grafting (CABG) and percutaneous coronary intervention (PCI)), was obtained from medical files. ResultsIndividuals with CAD exhibited higher total, periventricular, and deep WMH volumes than controls. Participants who underwent CABG had higher superficial WMH volumes than those with PCI, suggesting greater disease severity influences WMH burden. ConclusionCAD is characterized by a distinct pattern of cerebrovascular vulnerability, with revascularization procedures influencing WMH burden

Spousal caregivers of individuals with Alzheimers disease and related dementias frequently experience elevated perceived stress, caregiver burden, and loneliness, which are associated with adverse health outcomes. Early identification is therefore critical for timely intervention. Existing approaches commonly rely on wearable sensor data and standardized psychological questionnaires, while recent multimodal methods aim to improve prediction by integrating behavioral and linguistic information. In this study, we explored three modality configurations, wearable-derived features, interview-based text, and their combination, to classify caregiver psychological risk using the Perceived Stress Scale (PSS), Zarit Burden Interview, and UCLA Loneliness Scale. We compared traditional machine learning models and large language models (LLMs) (Gemini 2.0, Llama 4, and GPT-4o) under psychometrician-centered and caregiver-centered prompting strategies. Traditional machine learning models performed better under multimodal settings, while LLMs achieved stronger performance with Interview-Only input. We further demonstrate that PSS was the most predictable construct and prompting strategies substantially influenced LLM performance.

Major depressive disorder (MDD) is a common psychiatric disorder that requires reliable and objective assessment for early clinical intervention. Electroencephalography (EEG) is widely used for this purpose because it provides a non-invasive and low-cost measure of brain activity with high temporal resolution. However, EEG-based depression detection remains challenging due to the nonlinear nature of EEG signals, inter-subject variability, and the limited availability of subject-independent evaluation. To address these issues, this paper proposes a hybrid quantum-classical multiscale long short-term memory with parameterized quantum circuit branches (MS-LSTM-PQC) framework for subject-level EEG-based depression detection. The proposed model extracts temporal representations at multiple scales using parallel LSTM branches and incorporates eyes-closed (EC) and eyes-open (EO) condition information through condition-aware feature fusion. To further enhance the learned representations, scale-specific LSTM features are processed using PQC-based quantum branches implemented with TensorFlow Quantum (TFQ), providing an additional nonlinear feature transformation before classification. Experiments were conducted on the Mumtaz EEG depression dataset using EC-only, EO-only, and merged EC+EO conditions with 1-s, 2-s, and 3-s EEG windows. To reduce subject-level data leakage, all experiments were evaluated using 5-fold and 10-fold GroupKFold validation. The best overall accuracies across the evaluated settings were 92.05% and 95.08% under 5-fold and 10-fold GroupKFold validation, respectively. The 2-s merged EC+EO setting provided the most stable performance across validation protocols. In addition, Integrated Gradients (IG)-based explainability analysis showed that frontal and fronto-central channels, especially Fz, showed higher contributions to the model decision. These results suggest that multiscale temporal learning with quantum-enhanced feature transformation can support subject-level EEG-based depression detection under leakage-controlled evaluation.

The choroid plexus (ChP) plays a central role in cerebrospinal fluid production, immune signaling, and metabolic clearance, and has emerged as a potential imaging biomarker of neurodegeneration. However, accurate and scalable quantification of ChP volume remains challenging due to its complex morphology and low contrast on conventional MRI. The Automatic Segmentation of Choroid Plexus (ASCHOPLEX), a deep learning framework originally trained on healthy controls and multiple sclerosis cohorts, has not been systematically evaluated in neurodegenerative populations. Using T1-weighted MRI from the multi-center COMPASS-ND study, we assessed standard ASCHOPLEX performance in cognitively unimpaired (CU), Alzheimer's disease (AD), and Parkinson's disease (PD) participants (N = 30), followed by fine-tuning using expert manual segmentations (N = 60). Segmentation accuracy was evaluated using Dice, Jaccard, precision, and recall. The fine-tuned model was then applied to a larger cohort (N = 277) to derive normalized ChP volumes, which were compared across diagnostic groups using linear regression models. Fine-tuning significantly improved segmentation accuracy across all metrics (Dice: 0.45 to 0.84; Jaccard: 0.32 to 0.73; all p < 0.0001), enabling robust ChP delineation across sites and conditions. In the full cohort, normalized ChP volume was significantly higher in AD compared with CU and PD (p < 0.0001), while PD did not differ from CU (p = 0.31). These findings demonstrate that dataset-specific adaptation is essential for deploying deep learning segmentation models in heterogeneous neuroimaging cohorts. The refined ASCHOPLEX framework enables scalable ChP quantification and supports its use as a structural imaging marker in neurodegenerative disease.

BackgroundWhile resting-state fMRI demonstrated that brain networks are spatially dynamic (expanding, shrinking, and changing complexity over time), understanding the transient spatial network interactions that remain poorly characterized is critical for revealing the mechanisms underlying brain disorders. MethodsWe introduce DESINE (Dynamic Estimation of Spatially Interactive Networks), a novel framework using joint density distributions (2D histograms) of voxel-wise activity to quantify 4D spatial network interactions across sliding windows. We analysed transient deviations from the average functional state using root-mean-square error (RMSE) and mean absolute deviation (MAD), and characterized recurring interaction patterns using k-means clustering. We applied DESINE to 91 network pairs (14 networks) in a cohort of 508 subjects (315 healthy controls; 193 patients with schizophrenia, SZ). ResultsSZ is characterized by a significantly "constrained dynamic repertoire" of network interactions. SZ patients showed markedly lower means and standard deviations for both RMSE and MAD metrics across network pairs, particularly in regions of high activity, indicating systematic rigidity. Cluster analysis revealed significant alterations in state affinity metrics, suggesting a global breakdown in the brains capacity to preserve diverse, high-fidelity spatial configurations. Critically, these interaction metrics were associated with cognitive performance, symptom scores on the positive and negative syndrome scale, and chlorpromazine equivalent drug scores. ConclusionsThis work introduces DESINE as a global, voxel-agnostic framework for characterizing time-varying spatial interactions. Our findings highlight spatial rigidity as a fundamental feature of psychopathology, suggesting that the inability to express a diverse range of spatial interactions is a factor underlying cognitive deficits in schizophrenia.

Diffusion-weighted magnetic resonance imaging (dMRI)-based tractometry enables the quantification of white matter tissue properties in living humans while preserving anatomical specificity. Although tractometry is highly reproducible when the same scanner and acquisition protocol are used, its generalizability across scanners and protocols remains unclear. To address this gap, we performed a traveling-head experiment involving five subjects to evaluate tractometry across progressively different acquisition conditions, including multiple scanners, different scanner models, and two distinct protocols. Tractometry was performed for 20 major white matter tracts using diffusion tensor imaging metrics, neurite orientation dispersion and density imaging (NODDI) metrics, and a semi-quantitative ratio metric (T1w/b0). Generalizability across dataset pairs was quantified using the intraclass correlation coefficient (ICC). Tractometry showed consistently high ICCs when the scanner and protocol were identical; however, ICCs declined as differences in scanner model and acquisition protocol increased. Fractional anisotropy and orientation dispersion index retained relatively high ICCs across these comparisons, whereas other metrics showed marked declines when scanners or protocols differed. ComBat harmonization partially mitigated these declines, but ICCs did not reach the levels observed for datasets acquired using identical scanners and protocols. Finally, the minimum detectable change (MDC) for tractometry in datasets pooled across scanners and protocols varied by tract; for example, the optic radiation showed a lower MDC than the cingulum hippocampus. These findings highlight both the strengths and limitations of tractometry in multisite studies and highlight the importance of quantifying scanner- and protocol-dependent effects for specific metrics and tracts when interpreting measurements from heterogeneous datasets.

Brain similarity networks (BSNs), extracted from structural magnetic resonance imaging, provide a validated framework for studying brain network organization and encode neurodevelopmental information relevant for psychiatric disorders. Recently, a neurodevelopmental hypothesis has been proposed for bipolar disorder (BD), where evidence demonstrates neuroprogression phenotypes differing from controls. BSNs offer a promising framework for investigating BD's neural correlates but remain largely underexplored. Parallelly, graph neural networks (GNNs) have emerged as suitable deep learning models for exploiting network-level information. This study aimed to investigate BSNs for discriminating subjects with BD from controls within a GNN framework using the multi-site StratiBip network, composed of 605 controls and 501 subjects with BD. Leveraging advanced analysis tools, we developed a multi-site classification framework including: i) the state-of-the-art MIND algorithm for computing morphometric similarity (MS) networks based on gray matter volumes (GMV), ii) MS integration with age, sex, and GMV, iii) a leave-one-site-out cross-validation for multi-site model generalizability evaluation. The best model achieved a mean multi-site accuracy of 68%. Explainability analyses revealed meaningful MS patterns in the basal ganglia, frontal and temporal lobes, and a particularly relevant integration with age. This study provides interpretable insights into the role of MS in BD and unveils evidence supporting ageing-related processes as a significant component of BD pathophysiology.

Spontaneous thought is pervasive in everyday human cognition, yet datasets capturing its neural dynamics under minimally interrupted conditions remain limited. The current dataset was acquired from a think-aloud functional MRI experiment in which 118 participants continuously verbalized their spontaneous thoughts during 10-minute scanning sessions. The raw MRI data and verbal transcripts with sentence-level timestamps were previously released and analyzed in our prior study examining neural activity associated with thought transitions. Building on that release, we additionally provide preprocessed MRI data, speech transcriptions with word-level timestamps aligned to image acquisition, large language model-generated ratings of transcribed thoughts across emotional and sensory dimensions, and self-report survey measures assessing personality, mental health, and cognitive abilities. Validation analyses demonstrated activation in expected cortical regions associated with speech production and sensory content identified from transcript annotations, agreement between language model and human ratings, and adequate internal consistency of survey measures, supporting the datasets overall quality. This dataset enables reuse for investigations of spontaneous thought, speech generation, and individual differences using naturalistic functional MRI data.

IntroductionWhile global topological properties of brain networks reach relative maturity early in development, functional reconfigurations at the regional level continue throughout adolescence to support cognitive maturation. However, regional age and sex-specific developmental patterns of functional reconfiguration remain incompletely understood. MethodsWe analyzed resting-state fMRI data from 528 participants aged 5-21 years from the Human Connectome Project in Development. Three regional graph-theory metrics (betweenness centrality, hub score, and local efficiency) were computed for each individuals functional network. Cognition was measured using NIH toolbox. Parallel factor analysis was employed to decompose an individual x region x metric array into factors representing distinct developmental properties in the full sample and separately for males and females. Brain-cognition associations were examined in developmental subgroups (<13, 13-18, >18 years). ResultsThree factors emerged, characterizing visual, multimodal integration, and higher-order factors. Across development, metrics capturing network integration (betweenness centrality and hubness) showed general stability, while metrics capturing segregation (local efficiency) presented distinct peaks, particularly in the visual factor. Females showed earlier peaks and declines in higher-order factor, while males exhibited greater variability and protracted maturation in multimodal and higher-order factors. Brain-cognition associations were modest with early childhood and crystallized cognition composites showed small negative correlations with hub score in entire sample (r=-0.212) and local efficiency in males aged <13 years (r=-0.215). ConclusionFindings highlight nonlinear, sex-specific functional reconfiguration at region-level during childhood and adolescence, underscoring the importance of sex-stratified analyses in developmental and providing a crucial foundation for future investigations of developmental disorders.

Generalizable neuroimaging biomarkers that detect cerebral cortical changes after traumatic brain injury (TBI) and predict patient outcomes are needed to improve care and to develop targeted therapies. We used morphometric inverse divergence (MIND) analysis of structural MRI to investigate cortical gray matter morphological networks cross-sectionally and longitudinally after TBI and correlate these with symptoms, disability and cognition six months after injury. Our findings support the Triple Network Model from functional MRI of post-traumatic alterations in the relationship between task-positive, default mode and salience networks. However, the strongest associations between early cortical similarity metrics and long-term patient outcomes involved the dorsal attention network and the limbic network as well as similarity metrics across Mesulam's hierarchy of laminar differentiation. Since MIND mapping of cortical gray matter networks only requires data that is a routine part of standard clinical MRI protocols and does not need image harmonization across different scanners, this work reports a promising new tool that is immediately available for advancing research and clinical care in TBI.

Hypoxia-targeted BOLD MRI is a novel technique, which probes oxygenation physiology in response to a controlled transient hypoxia stimulus. In glioblastoma, the signal response is spatially and temporally heterogeneous. We developed a voxel-wise temporal decomposition framework for hypoxia-targeted BOLD MRI that separates the arrival of responses, transition phases, and steady state during controlled isocapnic hypoxia. Twenty healthy controls underwent 3-T BOLD MRI during a double hypoxic step challenge to establish a normative reference. Three patients with newly diagnosed glioblastoma were included as proof-of-concept cases. For each voxel, we estimated response arrival delay (Delaycorr), delay to plateau, delay to return and an O2-normalized steady-state response (HypoxiaSS). Healthy-control maps were used to construct a voxel-wise normative atlas and, for HypoxiaSS, a global-response-adjusted model for patient deviation mapping. In healthy controls, HypoxiaSS showed lower supratentorial between-subject variabilitythan both whole-stimulus comparators (coefficient of variation: 1.77 versus 2.36 for Hypoxiaavg) and higher voxel-level step-to-step agreement (ICC(2,1): median 0.951 versus 0.792 for Hypoxiaavg). Whole-stimulus averaging exhibited a systematic step-2 signal amplification present in 19 of 20 subjects, which was absent from HypoxiaSS. Asingle global response scalar explained a median 72.5% of voxel-wise between-subject variance in HypoxiaSS. In proof-of-concept patient analyses, G-adjusted HypoxiaSS deviation maps and timing maps identified spatially coherentabnormalities that were partly complementary and extended beyond conventional MRI-defined lesion margins.Temporal decomposition improves the stability and interpretability of hypoxia-targeted BOLD MRI and provides a practical framework for population-referenced physiological mapping and atlas-based deviation mapping in glioblastoma.